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Abstract. We show that the geodesic diameter of a polygonal domain 
with n vertices can be computed in 0{n* log n) time by considering O(n^) 
candidate diameter endpoints; the endpoints are a subset of vertices of 
the overlay of shortest path maps from vertices of the domain. 

1 Introduction 

For many geometric shortest path problems efficient solutions have been devel- 
oped both for simple polygons and for polygonal domains with holes, as well 
as for surfaces of polytopcs in 3D [9]. A notable exception is computing the 
diameter. The problem is non-trivial already in simple polygons, where it was 
examined decades ago [4,12] culminating in a linear-time algorithm [7]. Simi- 
larly, for convex polytopcs polynomial-time algorithms have been known since 
1990's [1,10]; the current best running time is 0(ri^ logn) (for a polytope with n 
vertices) [6]. However, for polygonal domains with holes, no algorithms existed 
until very recently. (The survey [9] mentions "brutc-forcc" results, but no details 
are given.) 

In [3], Bae et al. presented an 0{rJ^''^^^ logn)-timc algorithm for computing 
the diameter of an n- vertex polygonal domain P.^ Each of the diameter end- 
points may be a vertex of P, a point in the relative interior of P's edge, or an 
interior point of P. If cither of the endpoints is a vertex, the diameter can be 
found in O(n^logn) by computing the shortest path maps [8,9] from all ver- 
tices. However, Bae et al. [3] exhibit examples where the diameter is realized by 
non-vertex points on the boundary of P or by a pair of interior points. Accord- 
ing to [3], "This observation also shows an immediate difficulty in devising any 
exhaustive algorithm since the search space like dP or the whole domain P is 
not discrete." 

In this paper, we show that there actually docs exist a discrete set of can- 
didate diameter endpoints: an 0(n^)-size subset of vertices of the overlay of 
shortest path maps from vertices of the domain. For each candidate, the far- 
thest point of P can be found in 0(n log n) time by building the shortest path 
map. 

^ Earlier, in [2], the same authors announced an 0(n®"'"TT+^)-time algorithm. 



2 Preliminaries 



Let P be a polygonal domain with n vertices. The shortest path map from a 
source point p E P [9] is the decomposition of P into cells with the same unique 
combinatorial type of the shortest paths from p to the points in a cell (the 
combinatorial type of a path is the sequence of vertices visited by the path). 
We use SPM as a shorthand for " shortest path map" ; the SPM from a specific 
source p is denoted by SPM(j3). A bisector in SPM is the boundary between 
two cells, i.e., an edge of the SPM. The vertices of the SPM are of three types: 
vertices of P, intersections of bisectors with sides of P, and triple points where; 3 
or more bisectors meet. The complexity of the SPM is 0{n), and it can be built 
in O(nlogn) time [8]. 

Let A be the overlay of the SPMs from all vertices of P. The overlay can be 
built in 0(71^) time since there are 0{n) edges in the SPM from each of the n 
vertices, and two edges intersect 0(1) times. We will use the term node for a 
vertex of A, to distinguish nodes from vertices of P. 

Let d{p,q) denote the geodesic distance within P between points p,q, and 
let s*,t* be the endpoints of the diameter of P: d{s*,t*) = ma.Xp^q^pd{p,q). As 
in [3], we separately consider the cases when both of s* , t* are interior, when one 
of s*,t* is on the boundary, and when both are on the boundary. The hardest, 
bottleneck case is when both diameter endpoints are interior. In this case, the 
algorithm of [3] makes use of the following: 

Lemma 1. [3] If both s* and t* are in the interior of P, there are at least 

5 homotopically different shortest s* -t* paths. In addition, there exist at least 
3 distinct vertices ui,U2,U3 that are adjacent to s* on the shortest s*-t* paths; 
similarly, there are at least 3 distinct vertices vi,V2,V3 adjacent to t* on the 
shortest paths: d{s*,t*) = \s*Ui\ +d{ui,t*) = \t*Vj\ +d{vj,s*) for i,j = 1,2,3. 

The above properties are analogous to the properties of the diameters of polyhe- 
dral surfaces in M'^, established by O'Rourkc and Schevon [10]. For a formal proof 
of the properties refer to [3]. An intuitive explanation of the 5-path property is 
that each of s* , t* has 2 coordinates (degrees of freedom) , and the equation that 
two shortest paths are the same length takes away one degree of freedom; thus, 
5 path length equalities (i.e., 4 equations) take away all 4 degrees of freedom, 
pinning s* and t* . The 3- vertices property is even more intuitive: were there only 
2 vertices adjacent to s* through which the shortest s*-t* paths go, one could 
move s* away from both of them, increasing the diameter. 

Based on Lemma 1, Bae et al. [3] scroll through all 5-tuples of vertices of P, 
and for each 5-tuple look at each of the O(n^) cells in the overlay of the SPMs 
from the 5 vertices; within one cell, a constant-size system of constant-degree 
equations is solved to obtain candidate diameter endpoints. For each of the 
0{n'^) candidate pairs of points, a two-point shortest path query is performed 
in 0(n*/^^ log n) time [5]; the data structure for the queries can be built in 
0{n'^~^^^^^~^^) time which is the ultimate running time of the algorithm of [3]. 



Our algorithm identifies an 0(n^)-size set of candidate diameter endpoints. 
For each candidate, we simply find the farthest point in P with the O(nlogn)- 
time algorithm of [8]. 

2.1 Overview of our approach 

We say that a vertex u is the first bend of an s*-t* path if the path starts 
from the segment s*u; similarly v is the last bend if the path ends with vt* . 
(Recall that none of s*,t* is a vertex of P.) Wc start from proving a little 
variation of Lemma 1: If both s* ,t* are interior, then there exist exactly 3 vertices 
{wi, W2, Us}, exactly 3 vertices {v\,V2, v^}, and at least 5 homotopically different 
shortest s*-t* paths such that {'^1,^2,^3} are the first bends of the 5 paths, 
and {vi,V2,V3} are the last bends. (We only need the existence of the vertices 
and the paths; algorithmically we do not scroll through all triples of vertices as 
potential bend points.) 

We then define a graph G on s* ,t* ,u-i,U2,U3,Vi,V2, V3 that reflects the way s* 
and t* are connected with the 5 paths. A simple case analysis of the connectivity 
oiU = {ui, U2, U3} and V = {vi, V2, V3} in G shows that each of s*, t* is cither a 
triple point in the SPM from a vertex of P, or is at the intersection of bisectors 
in SPMs from two vertices (i.e., at a node of A). 

Moreover, in the latter case, the vertices of G are in a special relation; we 
prove that there are only 0{n^) sets of vertices that can be in the relation. Thus 
there are only O(n^) candidate diameter endpoints; for each, farthest point in 
P can be found in O(nlogn) time. 

2.2 Triple-point-diameter-end — an easy case 

Suppose that there exists a diameter whose (at least one) endpoint is a triple 
point in the SPM from one of the vertices of P; wc say that such point is a triple- 
point- diameter- end. In this case the diameter can be computed in O(n^logn) 
time since the total number of triple points in all SPMs in O(n^), and for each 
triple point the farthest point of P can be found in 0(n log n) time. In what 
follows we will assume that there is no diameter, one endpoint of which is a triple- 
point-diameter-end, since if it is the case we can compute it in 0{n^ log n log n) 
time. 

2.3 An imaginary perturbation 

Lemma 1 states that the set of the first bends of the (at least 5) shortest paths 
defining the diameter, has at least 3 vertices. The statement does not exclude the 
possibility that there exist, say, 5 distinct vertices each being the first bend of 
one of the 5 (or 6, or 7, or more) different shortest s*-t* paths. (In fact, it is easy 
to concoct instances with an arbitrarily large set of possible first bends of the 
diameters.) The next lemma shows that there actually exist exactly 3 vertices 
that serve as first bends for at least 5 shortest s*-t* paths. 



Lemma 2. There exist 3 vertices ui,U2,U3 and 5 homotopically different short- 
est s*-t* paths such that the set of the first bends of the 5 paths is {mi, U2, U3}. 

Proof. By Lemma 1, the set S of the first bends of the shortest s*-t* paths 
contains at least 3 vertices; clearly, s* is in the convex hull of S (or else s* can 
be moved increasing the diameter). Suppose that there are actually more than 3 
vertices in S. By Caratheodory's theorem, all but 3 points in S are "redundant" 
in the sense that there exist 3 points Ui, U2, U3 € S such that s* is in the convex 
hull of Ml, U2, W3. 

Let u be a vertex of P in 5 \ {ui, M2, U3}. Let P' be P with u replaced by 
a little notch so as to increase the length of paths that bend at u (Fig. 1); let 
q) be the geodesic distance within P'. Because, clearly, u is not on a shortest 
path from t* to any of ?ii,U2,U3, we have that (i'{s*,t*) =d{s*,t*). 

Claim. s*-t* is a diameter of P' . 

Proof. Suppose that there exists a pair of points a,b £ P' such that d'(a, 6) > 
d'{s*,t*). Then all a-b paths must go through u. (Otherwise, if an a-b path tt 
does not go through u then the lengths of tt in P and P' arc the same, and thus 
d(o, 6) = d'(a, 6) > d'{s*,t*) = di{s*,t*) contradicting the fact that s*-t* is a 
diameter of P.) If a-b were not a diameter in P (i.e., if d(a, 6) < d(s*, t*)), then 
we could perturb u by a small enough amount to make sure a-b is also not the 
diameter in P'; thus, a-b must be a diameter in P. But since there arc at least 
3 paths from a to b, there are also at least 3 paths from u to a, i.e., a is a triple 
point in SPM(u), or is a triple-point-diameter-end - a contradiction. □ 

In P', none of shortest s*-t* paths goes via u; thus the set 5 of the first bends 
of s*-t* diameters has decreased. We can continue this way, decreasing S until it 
has only 3 vertices. By Lemma 1, there still exist 5 s*-t* paths (in the perturbed 
P). But none of the 5 paths uses any of the perturbed vertices; hence the paths 
are the same in the perturbed P as they were in P - with {ui, U2, W3} as the set 
of the first bends. □ 



Fig. 1. Perturb P to remove shortest s*-t* paths not going through one of wi, W2, ws. 



We emphasize that the perturbation in the above proof is imaginary, not al- 
gorithmic (symbolic, random, or otherwise). Our algorithm does not use the 
existence of exactly 3 or more than 3 first bend vertices. We need Lemma 2 only 
to prove the correctness of our solution; the algorithm itself docs not have to 
find the 3 vertices, nor to perturb P to P', nor to do anything else according to 
the lemma. 




2.4 5 "independent" paths are necessary 

The 5-paths property ensures that whenever (s* , ) is a " diametricahy maximal" 
pair [3] (i.e., local motion of s*,t* cannot increase the geodesic distance between 
them), there exist 5 shortest s*-f* paths. However, the converse docs not hold 
automatically. That is, just mere existence of 5 shortest s*-t* paths does not 
make {s*,t*) diametrically maximal; it is also important that no subset of the 5 
paths could be obtained as a concatenation of a smaller number of (sub)paths. 

In particular, suppose that there exist 3 shortest s*-t* paths 7ri,7r2,7r3 having 
different first and last bends ui,U2,us,vi,V2,V3. Moreover, suppose that 7ri,7r2 
intersect other than at s* ,t* . Of course, two shortest s*-t* paths cannot properly 
cross. By "intersection" we mean that 7ri,7r2 partially overlap, sharing part of 
the way; i.e. that the paths have at least one common vertex v (Fig. 2). 




Fig. 2. Thick segments are the obstacles. The 5 shortest paths are tti = 

S*-Ul-V-Vl-t* ,TT2 — S* -U2-V-V2-t* , TT3 = s" -U3-V3-t* , TT4 = S* -Ul-V-V2-t* , TT5 = 

s*-U2-v-vi-t* . The figure is not to scale. 

In this case, there exist two more shortest s*-t* paths: 774 = s* -U\-v-V2-t* and 
7r5 = s* -U2-v-v\-t* . Nevertheless, even though there exist 5 shortest 8*-t* paths, 
s*-t* is not a diameter: the 4 paths length equality, |7ri| = |7r2| — |7r4| = \'k^\, 
are essentially only two equalities: |s*Wi| + d(wi,t>) = |s*U2| + d(w2, v), \t*Vi\ + 
d{v\,v) = \t*V2\ + d{v2,v). These equalities "eat up" only one degree of freedom 
from each of s* ,t* (putting s* on the U1-U2 bisector and t* on the V1-V2 bisector 
in SPM(i))). Equating the common length of 771,772,774,775 to the length of 773, 
takes away another degree of freedom from the pair {s* ,t*). Still, the pair retains 
one degree of freedom, and hence, s* ,t* can be simultaneously locally moved so 
as to increase the diameter. 

For a formalization of the above degree-of-freedom argument, one may look 
at the proof of Thm. 2 (Case (II)) in [3]. The end of p. 6 in [3] considers the 
case when the number of functions whose equality define {s* ,t*) is less than 5; 
or equivalently, when the number of equations for s* and t* is less than 4. This 
is exactly our case, as we have 3 equalities for the paths lengths. Bac ct al. prove 
that in this case the pair {s* ,t*) cannot give a local maximum of the geodesic 
distance. 



We summarize the above discussion in the following lemma: 

Lemma 3. Suppose that among the 5 shortest s*-t* paths there exist two such 
that the first bends of the paths are different, the last bends of the paths are 
different, and the paths intersect (overlap) other than at s*,t* . Then s*-t* is not 
a diameter. 

3 The diametric schema and an 0(n^ log n) algorithm 

We now give algorithms to compute the diameter in 0(n^ log n) time based on 
Lemmas 2 and 3. In the next section, we will carry out a more careful analysis 
of the bottleneck case of this algorithm to reduce the runtime to 0(n^ logn). 

3.1 Both s* and t* are interior 

We first consider the case when both s* , t* are interior points of P (the bottleneck 
case). By Lemma 2, there exist 5 shortest s*-t* paths, and triples of vertices 
U = {ui,U2, Us} and V = {vi,V2, v^} such that all 5 paths have one of U as the 
first bend and one of V as the last bend. Consider the graph G on s* ,t* , U, V, 
with edges between s* and each of U, between t* and V, and between u U 
and V G V whenever one of the 5 paths goes via u and via v (Fig. 3). We call G 
the diametric schema because it shows how s* and t* are interconnected by the 
diameters. We define the degree of it G i7 to be the number of its neighbors in 
V; the degree of v gV is the number of v's neighbors in U. 




Fig. 3. (a) If u €U has degree 3, t* is a triple point of SPM(w). (b) If ui and U2 have 
the same 2 neighbors, G is not planar, (c) The neighbors of U2 are V2,V3; w.l.o.g., us 
is connected to vs. 

We will now go through all possible interconnection patterns between U and 

V. First, suppose that a vertex u E U has degree 3 (Fig. 3a). In this case, there 
are 3 homotopically different shortest paths from u to t*: via vi, via V2 and via 
V3. That is, t* is a triple-point-diameter-end. 

Assume now that no vertex in [/, V has degree 3 in the schema. Since the 
total degree of the vertices in U is 5, two of the vertices, say, Ui,U2, are degree-2. 
Let vi , V2 be the neighbors of wi . 

First consider the case when is not a neighbor of U2. That is, the neighbors 
of U2 are vi,V2, and hence, the neighbor of ^3 is V3 (Fig. 3b). Then G is not 
planar: by contracting the edges s*U3,t*V3, we obtain K^^s (with s*,Vi,V2 in 



one part, and t*,ui,U2 in the other) as the schema's minor. However, the edges 
s*Ui, s*U2, s*U3,t*Vi,t*V2,t*V3 arc pairwise-non-crossing. Hence, there exist two 
shortest s*-t* paths such that the first bends of the paths are different, the last 
bends of the paths are different, and the paths intersect (overlap) other than at 
s* ,t* (Fig. 2). By Lemma 3, s*-t* is not a diameter in this case. 

In the remaining case, V3 is a neighbor of U2- That is, only one of Vi, V2 (say, 
V2) is a neighbor of 1*2; the other neighbor of U2 is W3 (Fig. 3c). This leaves two 
possibilities of connecting M3 to V: either to vi or to W3. Both possibilities result 
in the same (up to isomorphism) subgraph of the schema on U,V - a. 5-edge 
path through UUV. W.l.o.g. we will assume that this path is U1-M1-U2-U2-W3-U3, 
as in Fig. 3c. 

We are almost done: it follows from the schema that s* is on a bisector in 
SPM(t;2) and is on a bisector in SPM(w3); i.e., s* is a node of the overlay A. 
Thus, we can go through all nodes of A, and find the furthest point of P for 
each. Since the overlay contains O(n^) nodes, we find the diameter in 0{n^ logn) 
time. 

3.2 Both s* and t* are boundary 

When both s*,t* are on the boundary of P, each of them has one degree of 
freedom: [3, Theorem 2] proves that there must be at least 2 vertices that serve 
as first bends on the shortest s*-t* paths, and 2 vertices that serve as last bends; 
moreover, to pin both s* and t*, there must exist at least 3 homotopically differ- 
ent shortest s*-t* paths. Similarly to Lemma 2, we can assume that there exist 
exactly two vertices wi, U2, exactly two vertices vi,V2, and 3 shortest s*-t* paths 
such that U = {ui,W2} is the set of first bends of the paths and V = {v-i,V2} 
is the set of the last bends. Hence, w.l.o.g. the diametric schema looks as in 
Fig. 4a. This means that s* is a vertex in SPM(w2) as it lies on the intersection 
of a bisector in the map and an edge of P. Because there are O(n^) such vertices, 
the diameter can be found in 0{n^ logn) time. 



Fig. 4. (a) The schema for the case of both s*,t* being boundary: s* is a vertex of 
SPM(v2). (b) s* is boundary, t* interior: s* is a vertex of SPM(vi). 



3.3 s* is boundary, t* is interior 




Finally, if s* is on the boundary of P and t* is in the interior, there exist two first- 
bend vertices ui,M2, three last-bend vertices v-i,V2,vs, and at least 4 shortest 



s*-t* paths - this is proved in [3, Theorem 2] and can also be seen by the degrees- 
of- freedom argument. Thus, one of wi, f2, U3 (say, vi) has degree 2 in the schema 
(Fig. 4b), and s* is a vertex in SPM(i;i) (foot of a bisector touching an edge of 
P). Because overall there are O(n^) vertices in the SPMs from vertices of P, the 
diameter can be found in O(n^logn) time. 

4 Plausible vertices and an 0(n^ log n) algorithm 

The bottleneck case in the algorithm given in the previous section is when both 

s* , t* are interior and the diametric schema is as in Fig. 3c. The running time 
turned 0(n^ log n) because we scrolled through all O(n^) nodes of A. However, 
according to the schema, s* cannot be at an arbitrary node: it is defined by 
the intersection of a bisector between ui,W2 in SPM(u2) and a bisector between 
U2, U3 in SPM(w3), while W2,W3 define a bisector in SPM(u2). We now show that 
there are only O(n^) nodes with such properties. 

Say that vertices i,j are neighbors in SPM(fc) if there is a bisector between 
i,j in'sPM(fc). 

Definition 1. A node p of A is plausible if there exist vertices ui,U2,U3,V2,V3 
of P such that 

— p is at the intersection of the bisector between ui,U2 in SPM{v2) and the 

bisector between U2,u^ in SPM{v-:i); 

— 1^2,^3 are neighbors in SPM(u2). 

Definition 2. A 5-tuple (mi, W2, ""3, '^'2, i^s) of vertices of P is plausible if 

— 7/1.7/2 are neighbors in SPM(v2): 

— «2, it3 are neighbors in SPM{v3); 

— V2,vs are neighbors in SPM{u2)- 

The number of plausible nodes is not larger than the number of plausible 5- 
tuplcs (because even if a 5-tuple is plausible, the bisectors between ui , U2 and 

U2,U3 may not intersect at all). 

Lemma 4. There are O(n^) plausible 5-tuples. 

Proof. Define an n x n x n "bisectors adjacency" array B = {hijk} as follows: 



Here and throughout the indices run from 1 to n. 

A 5-tuple (ui,W2,W3,W2,W3) is plausible iff 6„i„2„2 = feuguaus = ^vsv^u^ = 1- 
The number of plausible 5-tuples is 




1 if there is a bisector between vertices i and j in SPM(A;) 
otherwise . 
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Since each SPM has hnear complexity, 

^braij=0{n) for all j. 

ml 

That is, for any j, there are 0{n) pairs {mUD such that b^j^j = 1. Hence, the 
number of plausible 5-tuples is 

0(n) 0(n) 
3 ml ik j t=l ik t=l ijk 

^ E m^Jii) il2K,mi} = 0{n) ■ 0{n) ■ 0{n) = 0{n') . 

t=l ij kj 

□ 

Returning to our algorithm, for every node of A, we can test in constant time 
whether it is plausible by checking the corresponding entries in B (clearly, B 
itself can be filled as the SPMs from the vertices of P are built). For each of the 
0{n^) plausible nodes, we find the farthest point in 0(n log n) time, and hence 

we have: 

Theorem 1. The diameter of P can he found in 0(n^ log n) time. 
5 Conclusion 

We showed how to compute the diameter of a polygonal domain in 0(n^ logn) 
time. A faster algorithm for the problem would have to use new insights: in our 
algorithm, already computing the arrangement A takes 0{ri^) time. 

An interesting open problem is whether our ideas can be applied to diameters 
of polytopcs in M"^ [1,6,10]. Shortest paths on polyhedral surfaces do not bend at 
vertices [11]; the combinatorial type of a shortest path is the sequence of edges 
that it visits (the path is uniquely defined by the sequence due to the unfolding 
property - the shortest path becomes a line segment if the polytopc is unfolded 
along the edges in the sequence). The 5-diameters property holds for polytopes 
as well [10]. If s* belongs to the interior of a face / of the polytope, then the 
diameters bend on at least 3 edges bounding /. By a perturbation argument 
as in Lemma 2, there exists exactly three edges U\,U2,Uz of / at which the 5 
paths bend; similarly, there exist 3 edges, vi,V2,vz, of f*'s face at which the 5 
diameters bend. Hence, one can define the diametric schema. However, here the 
analogy between polygonal domains and polytopes seems to end: were are not 
aware of a notion of a SPM from a polytope edge. 

The existence of several homotopically different paths between the diameter 
endpoints suggests to study properties of the 2nd, 3rd, and in general, Kih 
homotopically different shortest paths between two points in a polygonal domain. 
Even though algorithmically we do not use the paths to compute the diameter, 
it seems interesting to study their combinatorial properties. How can the "Kth 
SPM" be represented and what is its complexity? 
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